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This talk presents some results relevant for lattice QCD at higher order in ChPT First we discuss 
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volume formula for it. The latter allows for an alternative determination of meson sigma terms. 
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chiral perturbation theory of masses and decay constants of the charged mesons. We present 
results for all relevant quark mass combinations for the cases with two and three sea quarks. 
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Figure 2: The finite volume corrections R u to (uu). 



1. Introduction 

In this talk we present some work useful for the extrapolation of lattice QCD numerical results 
to the physical volume and quark masses. At low energies, there exists an effective field theory 
approximation to QCD known as Chiral Perturbation Theory (ChPT) [jl|]. We can thus use this the- 
ory to fit to lattice QCD numerical results and then extrapolate to infinite volume and the physical 
quark masses. ChPT is based on the correct inclusion of the nonanalytic structure in all amplitudes 
due to the Goldstone bosons of spontaneous chiral symmetry breaking and all relations that follow 
from the chiral Ward identities. This is done in a systematic expansion in momenta, energies and 
quark masses. There exists an enormous body of work in ChPT. Meson masses, decay constants 
and vacuum expectation values at infinite volume are known to two-loop order, see ^ [|, f|] and 
references therein. In this talk we present the extension of these results in two directions, finite 
volume for vacuum expectation values and partial quenching for the masses and decay constants. 

2. (qq) at Finite Volume 

The work presented in this section will be published in At finite volume, particles can 
propagate around the world. Taking this into account in ChPT was first done in [^]. Many more 
quantities have since been calculated at the one-loop level. Here we would like to address this 
problem at the two-loop level. A major problem is to evaluate the necessary loop integrals at two- 
loop order which has not been done so far 1 . An alternative approach to finite volume was pioneered 
by Liischer [f7|]. Here the leading effect of a finite volume is taken into account by an integral over 
another amplitude. An overview of earlier finite volume work can be found in [@]. 

We would like to check how well Liischer's method works when the finite volume quantities 
are fully calculated. At one loop in ChPT the full results for many quantities can be derived using 
Liischer's method. We therefore want to study a quantity at two-loop level. 

There are a few quantities where the two-loop calculation does not involve irreducible two- 
loop integrals. One such example is precisely the value of (qq) as can be seen from the diagrams 

'Some partial results were presented by C. Haefeli at this conference. 
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shown in Fig. []]. (qq) at two-loops at infinite volume is known [^]. We have redone this calculation 
checking explicitly that the 0{d — 4) parts of the integrals do not contribute to the final results. 
The known one-loop integrals at finite volume are therefore sufficient. Numerical results for R q = 
((qq)v — °°) / (?<?)<*> ^ shown in Fig. || for the case of the up quark condensate. We have here 
expressed the masses appearing in the integrals as the lowest order masses, leading to very small 
&(p b ) corrections. 

The analogue of Liischer's mass formula is derived in [Qj. The two types of contributions are 
shown in Fig. [|. Those of Fig. ^a do not contribute in ChPT to parity-even operators like O = qq. 
Keeping also the non-leading contributions of the same propagator, see |§|], we obtain 



(0) v -(0). 



1 



o 16712 ^ Jml + ql 



(4>\o\<$>) 



16^^) 



(2.1) 



with = \[km§L for the contribution of a neutral scalar (j> to the finite volume dependence of 
(O). m(k) is the number of times k = Y,i=\3 n j appears in the sum over n. 

Note that in this case there is a universal factor which contains the matrix element (<p\0\<j)) 
rather than an integral as in the finite volume mass formula. This also means that the finite volume 
dependence can be used to obtain this matrix element. Unfortunately, the extremely small €?(p 6 ) 
corrections allow no meaningful test of the Luscher formula versus the full 0(p 6 ) calculation. 



3. Partially Quenched ChPT at Two Loops 

The extension of ChPT to the quenched and partially quenched case has been done at one loop 
by Morel, Bernard, Golterman and Sharpe [K]] and discussed in detail in [11]. We use here the 
version without the super-singlet <J>q. Basically one adds bosonic ghost quarks to cancel the effect 
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of valence loops and a separate set of sea quarks to the Lagrangian of QCD, shown schematically 
in Fig. |]. As a consequence the Lagrangian gets (approximately, see [jll]]) a graded symmetry 
SU (n v + n s \n v ) x SU (n v + n s \n v ), with n v the number of valence and n s the number of sea quarks, 
rather than the usual chiral symmetry. We thus expect the low-energy effective theory to be Chiral 
Perturbation Theory with this symmetry group instead. Most of the classification work and the 
infinity structure at two loops Jl2| ] can be carried over to this sector by simply replacing traces by 
super-traces. Exceptions to this are baryons and the Cayley-Hamilton relations. 

One important note is that since QCD is a continuous limit of partially quenched QCD (PQQCD) 
the low-energy constants (LECs) from ChPT can be derived directly from those of PQChPT, e.g. 

L\=Lf pq) /2 + Lf pq) . (3.1) 

The superscript (3pq) means the constants relevant for 3 sea quark flavour PQChPT. 

The calculations published so far are: The mass at two loops in PQChPT for the case of equal 
valence and equal sea-quark masses for three sea flavours §T$j\, the decay constants for three sea 
flavours [14] and the masses and decay constants for two sea flavours [|TJ]. The remaining mass 



combinations for the masses are in the process of being written up [|TJ]. That article will also 
contain a longer introduction to PQChPT at two-loop order. 

The expressions are extremely long compared to those of the normal case [fj]. The main 
reasons are: First, there are many more masses even compared to the isospin breaking case [|]]. 
This allows for many more mass combinations to show up. Second and most important, the neutral 
or diagonal sector propagators contain both single and double poles with complicated coefficients. 
These coefficients are ratios of differences of lowest order meson masses to various powers. These 
rational functions of quark masses satisfy many nontrivial relations which allowed to compress the 
formulas by significantly more than an order of magnitude, and bring them to a publishable size, 
but they still leave an added layer of complexity. 

The double poles in the propagators appear because PQQCD is not a well defined field theory 
but only exists as a thermodynamical model. This leads to the well-known quenched logarithms [Q]. 

In the remainder we present some numerical results at two-loop order for masses and decay 
constants. We present them as a function of quark masses, m,-, via the lowest order meson masses 

Xi = mf ) =2B omj . (3.2) 

The border of validity for ChPT is for the %i of order 0.25 ~ 0.30 GeV 2 . 

One can of course show any amount of figures for the various mass cases for decays and 
masses. We restrict to a few simple cases. First the 1+1 case, i.e. valence masses equah^i = %2 = 
X3 and sea masses equal: X4 = X5 = X6- We show results along the lines X4 = tan^^i and X4 
constant, shown in Fig. || in the valence/sea quark plane. The result for the mass at one-loop order 
is shown in Fig. |6[ Notice the presence of the quenched chiral logarithm along the line labeled A. 
The same result at two-loop order is shown in Fig. [7|. The corrections are much larger than at 
one-loop order. This is also the case for normal ChPT There are choices of the LECs possible 
where this is not so [[J, 17]. We also show results for the 1+1 case of the decay constant Jl4|]. 



Plots for the other mass cases can be found in the published papers. The analytical form of the 



expressions can be downloaded from [18] and numerical programs are available from the authors. 
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In conclusion for this part, we have obtained the two-loop formulae for charged meson masses 
and decay constants in PQChPT and shown that the corrections can be sizable. We hope that this 
work will be useful for the numerical lattice QCD community. 
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Figure 7: The relative correction to the masses at Figure 8: The relative correction to the decay 
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